Preprint typeset in JHEP style - HYPER VERSION 



CP3-07-17 



MadGraph/MadEvent v4: The new web generation 



Johan Alwall 

Stanford Linear Accelerator Center, Stanford University, 
Stanford, CA 94309, USA 
E-mail: 



alwallOslac . Stanford . edi; 



Pavel Demin, Simon de Visscher, Rikkert Frederix, Michel Herquet, Fabio Maltoni 

Centre for Particle Physics and Phenomenology ( CPS), 
Universite Catholique de Louvain, 

Chemin du Cyclotron 2, 1348 Louvain-la-Neuve, Belgium 



E-mails: simon . de . visscher@f ymi . ucl . ac . be . pavel . deminOf ynu . ucl . ac . be 



frederix@fyma.ucl . ac . be , mherquet@f yma. ucl . ac .be , maltoniOf yma. ucl . ac .be 



Tilman Plehn 

School of Physics, University of Edinburgh, 

James Clerk Maxwell Building, The King's Buildings, Edinburgh EH9 3JZ, Scotland, UK 
E-mail: tilman.plehn@cern. cb| 



David L. Rainwater 

Department of Physics and Astronomy, University of Rochester, 
Rochester, NY 14627, USA 



E-mail: rain@pas.roch.ester.edu 



Tim Stelzer 

Department of Physics, University of Illinois at Urbana- Champaign, 
1110 West Green Street, Urbana, IL 61801, USA 



E-mail: tstelzer@uiuc.edu 



- 1 - 



Abstract: We present the latest developments of the MadGraph/MadEvent Monte Carlo 
event generator and several applications to hadron collider physics. In the current version 
events at the parton, hadron and detector level can be generated directly from a web 
interface, for arbitrary processes in the Standard Model and in several physics scenarios 
beyond it (HEFT, MSSM, 2HDM). The most important additions are: a new framework for 
implementing user-defined new physics models; a standalone running mode for creating and 
testing matrix elements; generation of events corresponding to different processes, such as 
signal(s) and backgrounds, in the same run; two platforms for data analysis, where events 
are accessible at the parton, hadron and detector level; and the generation of inclusive 
multi-jet samples by combining parton- level events with parton showers. To illustrate the 
new capabilities of the package some applications to hadron collider physics are presented: 

I. Higgs search in pp — > H — > W + W~: signal and backgrounds. 
II. Higgs CP properties: pp — > Hjj in the HEFT. 

III. Spin of a new resonance from lepton angular distributions. 

IV. Single-top and Higgs associated production in a generic 2HDM. 
V. Comparison of strong SUSY pair production at the SPS points. 

VI. Inclusive VF+jets matched samples: comparison with the Tevatron data. 
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1. Introduction 



Accurate simulation of both signal and background will play a key role in making discov- 
eries at the LHC. A well-known example is given by supersymmetric models, one of the 
most studied and cleanest templates for physics beyond the Standard Model. Character- 
istic signatures, from large rate multi-jet plus missing Et events coming from squark pair 
production, to extra 6-jets and r's from Higgs production, have in general large Standard 
Model backgrounds that need to be measured from the data and/or well described by 
Monte Carlo's. At the same time, the most important distinctive features of the signal, 
such as, for instance, mass distributions, kinematic edges and angular correlations will be 
exploited not only to improve the signal over background ratio but also to identify the 
quantum numbers (e.g., spin, color) of the intermediate heavy states which decay. 

The need for better simulation tools has spurred an intense activity over the last 
five years, that has resulted in several important advances in our ability to accurately 
simulate hard interactions. At the matrix element level, these include the development 
of general purpose event generators, such as MadGraph |ij and MadEvent J2|, Com- 
pHEP/CalcHEP |, §, SHERPA § and WHIZARD §], high efficiency multiparton gen- 
erators which go beyond the usual Feynman diagram techniques, such as ALPGEN j?j and 
HELAC |p, as well as Monte Carlo's that include NLO corrections, such as MCFM || and 
MC@NLO |0|. 

An accurate simulation of a hadronic collision requires a careful integration of the ma- 
trix element hard process with the full parton showering and hadronization infrastructure 



1 11 , |l^] . Here too significant advances have been made in the development of matching al- 



gorithms such that by Catani,Krauss, Kuhn and Webber (CKKW) [13, [lj, 15], by Mangano 



(MLM) |l6[ and by Lavesson and Lonnblad [17] and in their comparison [|1^, 19 1. A break- 



through was also achieved by Frixione, Webber and Nason [10, 2C] who showed how to 
correctly interface an NLO calculation with a parton shower to avoid double counting, re- 
leasing the first event generator at NLO, MC@NLO. Further developments and refinements 



in these directions are ongoing [EU 22]. 



While automatization has not yet been achieved at NLO, at tree- level we are now 
in a position where the matrix elements can be calculated by several different tools, and 
integrated in multiple ways into different parton-shower /hadronization codes to provide 
multiple checks of the accuracy of the simulation. 

Each of the tools mentioned above was developed with a unique approach optimized 
to the meet the authors intentions. This diversity of approaches, that include a wide range 
of overlap in physics reach is a key element to strengthening the programs and providing 
confidence in their results. It is important to understand the philosophy/intentions of the 
authors to understand and utilize their code. 

The single underlying principle in the development of MadGraph and MadEvent was 
to develop a tool that would maximize the amount of time the physicist could concentrate 
on physics, and minimize the obstacles between having an inspired idea, and being able 
to compare it to experimental data. This originated first with MadGraph and then with 
MadEvent. At a time when many phenomenologists were spending enormous amounts 
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of time and energy performing important tree level calculations, MadGraph/MadEvent 
was able to automate the whole process, from the calculation of matrix elements to the 
generation of unweighted parton-level events, allowing physicists to concentrate on other 
pressing issues. 

This current release of the MadGraph/MadEvent package is the natural extension of 
this project both in the theoretical, and the experimental directions. While the original 
MadGraph had the standard model "hard-coded" this version includes several new models 
(MSSM, 2HDM, HEFT,...), as well as the capability for user defined models. It also 
includes the ability to seamlessly pass the events through a full hadronic simulation so 
they can be subjected to a complete detector simulation. Utilizing the strengths of the 
web, a physicists can now literally go from an inspired concept to a full event simulation 
with detector reconstruction with just a few clicks of the mouse. In addition, in order to 
facilitate the communication between theorists and experimentalists, the complete process 
from diagram generation to event analysis can be recreated from a simple text file which 
collects all the input cards used in the various phases of the simulation. 

In this work we briefly illustrate the new features of the code, focusing more on the 
structure and on general aspects than on the technical details. In Sections 2 and 3 we 
give an overview of the package. In Section 4 we discuss the models for which a dedicated 
implementation exists (SM+HEFT, MSSM, 2HDM) as well as a framework where new 
models can be implemented starting from the SM. In Section 6 our approach to generate 
multi-jet samples through matching is described. Section 7 collects the descriptions of 
some of the tools available that allow the users to handle the parton-level events in the Les 
Houches format for different purposes: from analysis to further simulation, to the decay 
of unstable particles maintaining the leading spin-correlations. In Section 8 we provide 
several examples, from SM measurements to the search for SUSY in inclusive signatures, 
of studies that can be performed with our package. We draw our conclusions and present 
our current line of research and development in the final section. 

2. The path to event generation 

The new structure of the MadGraph/MadEvent package is shown in Fig. |l|. Each of the 
four steps (code creation, parton-level event, hadron-level event and reconstructed object 
generation) are driven by input cards provided by the user. All cards have a comment 
section which describes in detail their content and syntax. Here we just mention their 
functionalities and how the generation proceeds, as outlined in Fig. |l[ For sake of simplicity, 
we consider the procedure to be followed by a web user, who can control the full generation 
process on one of the MadGraph clusters, via a web interface. The steps undertaken by a 
user working locally on his/her own computer, are the same, but in this case the simulation 
is initiated by calling scripts from the command line. In this case, the exact syntax to be 
used can be found in the corresponding README files and in our on-line documentation. 

The first phase is the creation of the code specific to the process that the user re- 
quests. This can be done directly by filling a web form or by editing the process card, 
proc_card.dat, to identify the process or processes to be included in the code, together 
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Figure 1: Flow chart of the MadGraph/MadEvent event generation. Each step is controlled by a 
specific card. Details are given in the text. 



with the model and the maximum order of the couplings to appear in the tree-level di- 
agrams. Labels can be specified to identify multi-particles (such as jets) which defines 
summations over physical particles. Examples of process cards are given in the Section [7|. 
Once completed, the process card is read by MadGraph and a process-specific code (pack- 
aged in a gzipped tar file) is created and can either be downloaded in order to run the 
simulation locally, or be run directly on the cluster using a web interface. 

Two cards are needed to perform the event generation. The first one is the parameter 
card, param_card.dat, which contain the numerical values of the necessary parameters for 
a specific model. The parameter card has a format compliant with the SUSY Les Houches 
Accord and it is dependent on the physics model. At variance with the previous version, 
secondary parameters and widths are no longer computed in MadEvent, but instead passed 
through the parameter card also for the Standard Model. Such a card can be generated with 
an independent code, which we refer to as a "Model Calculator". There are calculators 
available on the web for the SM, MSSM and 2HDM models, which, starting from the 
parameters in the Lagrangian (primary parameters) calculate all secondary parameters 
(such as masses, widths, and auxiliary parameters) needed by MadGraph to perform the 
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cross section integration. We stress that the calculation of all such parameters should 
be done by a Model Calculator since they must be consistently performed at tree level. 
The second card to be completed is the run card, run_card.dat, where all the information 
regarding the event generation (number of events requested, random seed,. . . ), the collision 
(energy, beam type, pdf's, scales) and the acceptance of the detector (p™ m , r? m ax, • • • ) is 
passed. Both cards can be generated using a web form, or uploaded. Once the cards are 
completed, the event generation can start. 

As a result of the run, the user will obtain a data file with the unweighed events in the 
Les Houches format, a set of plots of kinematical variables, and a ROOT file. 

The parton-level events can then be passed to Pythia by specifying in the pythia_card . dat 
the desired options. Events at the hadron level are then passed to PGS, whose parameters 
are set in the pgs_card.dat. On the web the Pythia and PGS cards can be completed 
at the same time as the parameter and run cards or at a later time. In the first case the 
generation goes all the way through the reconstructed objects at the detector level, and 
the information (events and plots) of intermediate stages of the simulation is kept. 

The MadGraph/MadEvent package can be used at different levels, depending on the 
user's needs/expertise. We have identified three main categories: 

• A web user performs the full generation (from the process specification to the Pythia 
or PGS events) over the Web, downloading only the input cards and the results 
(events and plots). No code download or installation is needed as all phases are 
handled via a web interface and the generation is performed on one of our clusters. 
Each user manages a personal database where all the codes/runs/events are stored 
for future use. This is the simplest level of accessing MadGraph/MadEvent, and it is 
suitable to both theorists and experimentalists who might want a quick determination 
of a cross section or a fair amount of events to conduct preliminary studies using the 
available physics models. Code creation and all other tools on the web page are 
accessible to anybody upon web registration. 

For a web user it is enough to understand the logic and the steps of the code, Fig. [I], 
and to correctly manage the information which is passed through the input cards. 
This is explained in detail in our on-line documentation. 

• A local user creates the code specific to the process of interest over the web and then 
downloads and runs it locally on her/his personal computer or cluster. This allows 
somewhat more flexibility (for instance, special selection cuts can be implemented, as 
well as different PDF's, dynamic scale choices, and so on), as well as the possibility 
for the production of large samples of events. In this case, it is suggested to run in 
parallel over a cluster with the PBS [^] or CONDOR [^4|] batching systems. 

• A developer downloads the full MG/ME package to exploit all functionalities locally. 
This is certainly the preferred option not only for the users making detailed and/or 
sophisticated simulation studies, but also for developers of new tools or models in the 
MG/ME framework. We underline in passing that experiments or larger collabora- 
tions that have at their disposal a computer farm running a batching system, such as 
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Dir name 


Content 


Cards 
Source 
SubProcesses 
Events 
bin 
lib 
HTML 


steering cards 
proc. independent source files 
proc. dependent source files and dirs for the subprocesses 
LHE event, plots and ROOT files 
executables and scripts (csh and Perl) 
libraries and PDF data 
web pages 



Table 1: Organization of the process directory (Template). When the code specific to a user's 
process is generated over the web, it is packaged into MadEvent . tar . gz which has the structure 
above. 



PBS or CONDOR, might also install the web server software, which manages all the 
generation steps (code creation, event generation, interfacing to Pythia and PGS, as 
well as ROOT-based analysis) and user databases automatically. The corresponding 
web server software and tech support can be obtained from the authors. 

The installation of the MG/ME package is straightforward since no external libraries 
are needed. The interested reader should refer to the README file contained in the main 
directory, obtained after decompressing the MG_ME_V4.X.Y.tar .gz archive file. The main 
directory is organized in subdirectories whose contents we briefly describe below. 

• MadGraphll : The source code (Fortran) and executable (mg2) for the matrix-element 
generator. mg2 generates not only a code for the matrix element squared for a given 
process, but also all subprocesses contributing to it. When interfaced with MadEvent 
it also generates the mappings needed for the phase space integration, including 
the Single-Diagram-Enhanced expressions. mg2_sa is the "stand alone" version of 
MadGraph producing the code for the matrix element squared only. 

• Template : The general structure of a process dependent code, as shown in Table [I]. 
To create the code corresponding to a new process, a copy of the full directory should 
be made (at the main directory level) and the script . /bin/newprocess invoked from 
it. 

• Models : It contains the available models (each of them is a directory) including the 
user model template (usrmod) . The user developing a new model makes a copy of the 
usrmod directory at this level, giving a name to it that identifies the model (which 
can later be set in proc_card.dat). 



HELAS : The HELAS library ^§ source code. 



DECAY : The source code to decay unstable particles directly from parton-level LHE 
files. 
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3. The Standard Model implementation 
3.1 Standard Model interactions 

The Standard Model of particles and interactions, based on the SU(3) C x SU(2)l x U(1)y 
gauge symmetry has been available since the first versions of both MadGraph and more 
recently of MadEvent. There is, however, one important differences w.r.t. the previous 
version of the package, regarding how the couplings of the models are handled. As was 
already mentioned in the previous section, the task of computing from the parameters 
in the Lagrangian (primary parameters) all the secondary parameters (masses, widths 
and dependent parameters) needed by MadGraph is left to an external program, the SM 
Calculator. The output of the SM Calculator is a parameter card, param_card.dat, which 
contains the numerical values of the main couplings (primary and secondary) of a specific 
model. The parameter card has a format compliant with the SUSY Les Houches Accord. 

A simple example is given by the EW parameters that characterize the gauge SU{2) l x 
U(1)y interactions and its breaking: in the Standard Model there are five relevant param- 
eters, a em , Gf, sinOyy , mz, raw of which only three are independent at tree level. Various 
schemes differing by the choice of the parameters considered independent are used in the 
literature. In the SM Calculator, the default is to take GF,mz,mw as inputs and derive 
a em , sin 6\v, but other choices are available. As a result a consistent and unique set of values 
of the couplings appearing in the Feynman rules is derived and used for the computation 
of the amplitudes. 

Another sometimes important feature of our SM implementation, is the possibility of 
distinguishing between the kinematic mass (pole mass) for the quarks and that entering in 
the Yukawa coupling definition (MS mass). For the latter, the user can choose to evolve 
the mass to the scale corresponding to the Higgs mass, which leads to an improvement of 
the perturbative expansion. 

Finally, we mention that various versions of the Standard Model are actually available 
for specific studies. For example, in the "minimal SM" (sm) the CKM matrix is diagonal 
while in the sm_ckm model a mixing between the first and second generation is allowed 
(Cabibbo angle). Another example is the sm_nohiggs model where the Higgs has been 
eliminated and the EWSB sector behaves as a non-linear sigma-model. 



3.2 Higgs Effective Theory 

The Higgs effective field theory (heft) model is an 'extension' of the Standard Model (SM), 



where the Higgs boson couples directly to gluons (and photons) |26|, 27, 28]. In the SM 
these couplings are present through a heavy (top) quark loop. For a not too heavy Higgs 
(m/i < 2mt), it is a good approximation to take the mass of the heavy quark in the loop 
to infinity 1 . This results in effective couplings between gluons and Higgs bosons. 



^or this approximation to hold, not only should the Higgs mass be smaller than twice the top mass, 
also all other kinematic variables, such as the transverse momentum of the Higgs boson, should be smaller 
than 2m t I29I. 
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The effective vertices can be derived from the effective dimension five Lagrangian 

1 

4 



Ch = - 7 g h G%G a ^, (3.1) 



where G" = d^A® — d u A a ^ + gf abc A b ^Af,. The coupling constant is given by 

9h = (3.2) 

Due to the non-abelian nature of the SU(3)c color group the effective vertices do not only 
include two, but also three and four gluons coupling to the Higgs boson. Since MadGraph 
can work only with three- and four-point vertices, the four-gluon interactions in the heft 
model are obtained by rewriting the QCD four-gluon interaction in terms of three-point 



vertices with an extra non-propagating internal tensor particle, T |j(J, 31]. This trick can be 
easily understood by noting that the usual (text-book) form of the four-gluon interaction 
is the sum of three terms, whose color and Lorentz structure correspond to 2 — ► 2 diagrams 
where a color octet tensor is exchanged in the s,t,u channels. With the introduction of 
this extra particle, the four-gluon-Higgs vertices can be reduced to diagrams with at most 
four-point vertices. To get the standard diagrammatic visualization of four-gluon and 
four-gluon-Higgs vertices it is sufficient to contract the T particle lines to a single point. 

The gluon couplings to a pseudo-scalar Higgs are also implemented. The name of 
the pseudo-scalar Higgs in MadGraph is h3 (i.e., the same as in the 2HDM and MSSM 
models). The effective dimension five Lagrangian for the pseudo-scalar Higgs coupling to 
the gluons is 

C A = l -9AG%G a ^ A , (3.3) 

where (5° is the dual of G® , G" = \e^ vpa G a pa . The effective coupling constant qa is 
given by 

« = £■ M 

The pseudo-scalar Higgs has only effective couplings to two or three gluons. The four- 
gluon-pseudo-scalar Higgs vertex is absent due to the anti-symmetry of the epsilon tensor 
e tLVpa . If a mixed Higgs with no definite CP parity is needed, it sufficient to change the 
couplings of the Higgs to the gluons. First generate the process with the SM Higgs, then, 
after downloading the code, change the coupling in the . /Source/Model/couplings . f file. 
The coupling constant is defined as a two-dimensional object, where the first and second 
elements are the CP-even and CP-odd couplings of the Higgs to the gluons, respectively. 
The HELAS subroutines automatically use the correct kinematics for odd-, even- or mixed 
CP Higgs's coupling to the gluons. At present, the implementation allows production of 
only one Higgs-boson. The effective couplings of two Higgs bosons to gluons are available 
in HELAS, but not yet included in the HEFT model. 

4. Going beyond the Standard Model 

In the current version a few new physics models have been added and fully tested: the 
minimal supersymmetric extension of the Standard Model (MSSM) and the general two 
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Higgs doublet model (2HDM). In addition, a framework for setting up new models is 
available. Here, we briefly describe the main features of these new implementations. 



4.1 The 2HDM implementation 

The two-Higgs-doublet model (2HDM) has been extensively studied for more than twenty 
years, even though it has often been only considered as the scalar sector of larger models 
like the MSSM f3§ or Little Higgs models [§3||. The generic 2HDM considered here may 
display by itself an interesting phenomenology justifying its study. As a non exhaustive 
list, let us mention new sources of CP violation in scalar-scalars interactions [34|, tree-level 
flavor changing neutral currents (FCNCs) due to non diagonal Yukawa interactions, dark 
matter candidates [35] or Higgs bosons lighter than the LEP bound [36]. 

In the "full" version of the model (2hdm_f ull), no particular restrictions are imposed 
on the interactions allowed by gauge invariance, except electric charge conservation. Many 
diagrams involving tree-level FCNCs and violating the CP symmetry are thus present. 
The user who is not interested in these phenomena should use the "simplified" version of 
the model (2hdm), where the number of generated diagrams is in general much smaller. 

The following naming convention is used: h+ and h- stand for the positively and 
negatively charged Higgs bosons and hi, h2 and h3 stand for the neutral ones. Since the 
CP invariance of the potential is not assumed, the neutral bosons are not necessarily CP 
eigenstates and the standard naming convention in this case (i.e., hi being the lightest one 
and h3 the heaviest one) is used. 

TwoHiggsCalc is the calculator associated with the model. It has been written in C 
and is accessible from a web interface. It has been designed to compute input values for 
the 2HDM extension of MadGraph/MadEvent but it can also be used as an independent 
tool. Starting from various parameters of the Lagrangian, such as the vacuum expectation 
values (vevs) or the Yukawa couplings, the program computes useful secondary physical 
quantities at leading order such as the scalar mass spectrum, the mixing matrix, the total 
decay widths and the branching ratios. 

TwoHiggsCalc reads input and writes out results in a specific format close to the 
"SUSY Les Houches Accord 1.0" convention for SUSY parameters p7| . This format can 
later be read by MadEvent to perform numerical calculations for 2HDM processes. A 
README file describing this modified version of the LHA format used as input convention is 
available. To ease the use of TwoHiggsCalc, a web form has been designed to automatize 
the parameter card writing process. Numerical values for the parameters (units being 
fixed when needed) can be entered on this form. Some simple algebraic expressions can 
also be used. The +,-,*,/ operators and the reserved keyword PI, e.g., PI/2+3*PI/2, are 
correctly interpreted. 

In the general 2HDM, one has the freedom to choose a specific basis for entering pa- 
rameters. All the possible choices are physically equivalent (see e.g. |3^] for a discussion). 
TwoHiggsCalc and the 2HDM model both assume that the parameters are given in a partic- 
ular basis, called the "Higgs basis" where only one Higgs doublet gets a vacuum expectation 
value. An independent program , Gen2HB, has been written to convert parameters given 
in an arbitrary basis (where both Higgs doublets get vevs), called "generic", to parameters 
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in the Higgs basis. See [34] for more information on basis invariance and on the notation 
used. 

The scalar potential in the Higgs basis reads 

V = + V*hIh 2 - (mH\h 2 + h.c. N 



XiIhIhA + x 2 (h\h 2 



+A 3 (h\h 2 ) + A 4 (h\h 2 ) [Hl^ 

X 5 h\h 2 + XeHlHt + X 7 h\h 2 ) (h\h 2 ) + h.c. 



+ 

All parameters in front of quartic terms and the charged Higgs mass are input parameters, 
while fii, \i 2 and fi-j are fixed by minimization constraints and by the vev extracted from 
the observed SM parameters. Ai to A4 are real while A5 in general is complex. However, 
since only the phase differences between A5, Xq, A7 and /i 3 matter, the phase of A5 can 
always be rotated out. It is thus considered as a real parameter while Xq and A7 are a 
priori complex. 

In the same basis, the Yukawa interactions read 



r QlV2 
Ly = 



V 



{M d H x + Y d H 2 )d R + (M U H X + Y u H 2 )u R 



+ El^l [{ M e H 1 +Y e H 2 )e R ] . 

v 

Yukawa couplings are expected to be given in the physical basis for fermions, i.e., in the 
basis where the mass matrix is diagonal. Since in the Higgs basis only the first Higgs 
doublet gets a non zero vev, the M matrices are completely fixed by the physical fermion 
masses and CKM mixing matrix (restricted to Cabibbo angle) while the Y matrices (giving 
the couplings of the second Higgs doublet) are a priori free. For these matrices, the first 
index refers to doublet generation while the second refers to the singlet generation. For 
example, Y2B stands for the complex Yukawa couplings of the second Higgs doublet to the 
second generation quark left doublet and to the bottom singlet. 

In the generic basis, similar expressions are assumed. For the scalar potential all 
parameters in front of quartic terms are inputs as well as tan(/3), the norm of ^3 and the 
phase of v 2 . The overall vev is again extracted from SM parameters while mass terms 
parameters, like /Ui, [i 2 and the phase of // 3 , are fixed by the minimization constraints. Ai 
to A4 are real parameters, A5, Xq and A7 are a priori complex. Like in the Higgs basis, 
the Yukawa couplings must be given in the physical basis for fermions. Since the mass 
matrices are fixed, only the Yukawa coupling matrices of the second Higgs doublet (r), is 
required. The other one is going to be automatically evaluated to match observed fermion 
masses and CKM mixing matrix (restricted to Cabibbo angle). For the T matrix, the first 
index refers to doublet generation while the second one refer to the singlet generation. For 
example, G2B stands for the complex Yukawa couplings of the second Higgs doublet to the 
second generation quark left doublet and to the bottom singlet. 

Given the above parameters and some SM parameters, TwoHiggsCalc computes the 
following quantities 
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Block 


Comment 


SMINPUTS 


worn 1 to 4, blvl parameters, see the SM section tor more details 


MGSMPARAM 


Extra block with smVw and Myy, see the dM section tor more details 


MGYUKAWA 


Yukawa masses used m the Yukawa couplings evaluation 


MGCKM 


The full CKM matrix 




Basis choice, must be 1 (Higgs basis) for MadEvent ! 


MINPAR 


Scalar potential parameters in the Higgs basis 


YUKAWA2 


Yukawa couplings of the second Higgs doublet 


MASS 


All SM particles masses, plus the five new Higgs boson masses 


TMIX 


The scalar mixing matrix 


DECAY 


For all the Higgs bosons, top, and Z 



Table 2: LHA blocks used in the 2HDM implemention 



• Scalar particles mass spectrum 

• Normalized mixing matrix of neutral scalars (called T in |?4| ) 

• Decay widths for all scalars as well as for W and Z bosons and the top quark. All 
widths are evaluated at tree- level using the same couplings as in MadEvent. Below 
threshold formulas are included for the scalar decays into two vector bosons and the 
one loop driven scalar decay into two gluons is also computed. 

The LHA blocks and parameters used by MadEvent are given in Table All blocks 
in the table are provided by TwoHiggsCalc. Note that if parton density functions (PDFs) 
are used in the MadEvent run, the value for a s at Mz and the order of its running is given 
by the PDF. Otherwise a s {Mz) is given by block SMINPUTS, parameter 3, and the order of 
running is taken to be 2-loop. The scale where a s is evaluated can be fixed or evaluated 
on an event-by-event basis like in the SM. 



4.2 The MSSM implementation 

One of the most popular extensions of the Standard Model is TeV scale supersymmetry. 
Supersymmetry solves the problem of quadratically divergent corrections to the Higgs 
boson mass by the introduction of new bosonic particles having the same couplings as the 
Standard Model fermions, and new fermions having the same couplings as the Standard 
Model bosons, thus cancelling the loop contributions to the Higgs mass to all orders. The 
Minimal Super symmetric Standard Model, MSSM, represents the minimal particle content 
for a supersymmetric extension of the Standard Model together with the maximum coupling 
space allowed by so-called "soft supersymmetry breaking terms" in the effective low-energy 
Lagrangean. These are constructed not to introduce new divergencies in any couplings, 
and therefore maintain the cancellations of quadratically divergent corrections to the Higgs 



mass. For an introduction to supersymmetry and the MSSM, see e.g., Refs. 4C]. 



The implementation of the MSSM particles and vertices into MadGraph II was made 
in Ref. jll], E^], following the conventions of Refs. [43] and [44|. Specifically, it is restricted 
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to the minimal supersymmetric model conserving i?-parity, without CP-violating phases 
and with diagonal CKM and MNS matrices. Higgs Yukawa couplings as well as mixing 
between right- and left-handed sfermions are implemented only for the third generation. 
However, no specific super symmetry breaking scheme is assumed, so the spectrum and 
couplings of the supersymmetric particles can be produced with any spectrum generator 
regardless of the assumptions going into its calculations. The spectrum and couplings of 
the particles are read through SUSY Les Houches Accord files p7|]. 

In order to consistently calculate decay widths and the dependent parameters, a model 
calculator (see Section |2|) for the MSSM is available. MSSMCalc takes a SUSY Les Houches 
Accord (SLHA) file |[7| from any Spectrum generator as input, and produces a MadEvent 
readable file, param_card.dat, with the missing Standard Model parameters, as well as 
decay widths for all supersymmetric particles (calculated at leading order by Sdecay pq] ), 
the Higgs particles and the top, and Z particles. Care has been taken to ensure that 
the parameters used in the calculation of decay widths are as similar as possible to the 
parameters used in MadEvent, since the correct total decay widths are vital to get the 
correct tree- level cross-sections for processes involving decaying particles. 

In the default run mode, MSSMCalc uses the Standard Model parameters given in the 
SUSY Les Houches accord (a em , Gf and Mz) to calculate the parameters sinQy/ and M^y, 
which are stored in a MadEvent specific block MGSMPARAM in the resulting param_card.dat. 
The b quark pole mass is calculated from the MS mass at 2- loop order. Another option is 
to extract the Standard Model parameters (and the vacuum expectation value ratio tan (5) 
from the chargino and neutralino mixing matrices, in order to ensure unitarity of ino-ino 
scattering at high energy. In this mode, also the Yukawa masses of the third generation 
fermions are extracted from the third generation sfermion mixing matrices. For a thorough 



discussion of this option, see section II C of [41|. 



The strong coupling a s is calculated in MSSMCalc using 2-loop renormalisation group 
running in the MS scheme, at the scale specified in the GAUGE block statement. The value 
used for the strong coupling g in the decay width calculations is stored for comparison in 
the block GAUGE, parameter 3. Note however, that the value of a s used in MadEvent is 
given by the choice of parton distribution function and the scale chosen in the run. 

If there are blocks missing in the SLHA file which are necessary for MadEvent, MSSM- 
Calc will produce a param_card.dat file containing error messages. 

The SUSY Les Houches blocks and parameters used by MadEvent are given in Table |. 
All blocks in the table should be provided by the user (and are indeed provided by most 
MSSM spectrum generators), except for the MGSMPARAM and the DECAY blocks which are 
produced by the parameter calculator MSSMCalc. Note that if parton density functions 
(PDFs) are used in the MadEvent run, the value for a s at Mz and the order of its running 
is given by the PDF. Otherwise a s (Mz) is given by block SMINPUTS, parameter 3, and the 
order of running is taken to be 2-loop. The scale where a s is evaluated is however always 
given by the "scale" parameter in the run_card.dat. 

4.3 The User Model 

Beside the models already implemented in the MadGraph structure, it is also important to 
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Block 


Comment 


SMINPUTS 


Except tor 5, the b quark Mb mass 


» if /-I c* TiltTt ADAM 

MGSMPARAM 


hxtra block with sin% and Myy, written by MblbMOalc 


MASS 


Til" rj.iL 1 l 

Including 5, the b quark pole mass 


NMIX, UMIX, VMIX 




CTnDMTV CDflTMTV O T A T TM T V 

biUFnlA, brSU InlA, biAUniA 




ALPHA 




HMIX 


Only parameters 1 (fi) and 2 (tan (5) 


AU, AD, AE 


Only the third generation parameter 3 3 


YU, YD, YE 


Only the third generation parameter 3 3 


DECAY 


For all SUSY particles, Higgs bosons, top, W ± and Z 



Table 3: SHLA blocks used by SUSY MadEvent. See the text for details. 



let the possibility to the non-expert to be able to test his or her own models. For example, 
adding a Z' boson or a t fermion to Standard Model and all their possible couplings to 
other particles should be now straightforward. For this purpose a new framework has been 
developed to provide users with an easy and safe, yet flexible, tool for implementing their 
own models. It is important to mention that the user model is limited by the MadGraph 
and HELAS assumptions for the Lorentz and color structure of the vertices. 

The aim here is not to give a detailed step-by-step description on how to implement a 
new model (which is available in the README file in the usrmod directory), but to present 
the philosophy adopted in the design of the framework. In practice, a model is defined by 
a list of particles with given quantum numbers and their interactions. The implementation 
of a model proceeds in three simple steps. 

First, three files have to be provided: i) the list of particles; ii) the list of the interactions 
(the complete Standard Model, which is already present by default, plus the user's new 
inputs) iii) the list of possible new parameters present in non-SM couplings expressions. 
As a second step a Perl script, ConversionScript.pl, is invoked that takes the above 
information and creates all the files needed by MadGraph/MadEvent. As a final step, the 
user must modify by hand the couplings. f file to set the values of the couplings and 
provide the formulas to go from the parameters of the model to those appearing in the 
Feynman rules. 

Let us consider the implementation of a Z' vector boson as an example. In the particle 
list file (particle.dat), an extra Z appears as 

#Name anti_Name Spin Linetype Mass Width Color Label Model 
#xxx xxxx SFV WSDC str str STO str PDG code 

#MDDEL EXTENSION 

zp zp V W ZPMASS ZPWID S ZP 32 

where the needed information is provided in a simple syntax. In this example we have a 
Z' (zp) of vector type (V), with mass and width ZPMASS and ZPWIDTH. We made it a color 
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singlet by setting color to S. The linetype W and the label ZP are used in the Feynman 
diagrams for cosmetics only. 

As described above, in the second file (interactions.dat) the interactions are spec- 
ified. In this file the names of the particles as defined above have to be used. In this 
example the vertex Z'dd, i.e., the vertex between the Z' and two down quarks, is shown. 

# USRVertex 
d d ZP GZPD QED 

The name GZPD is choosen arbitrarily (but paying attention not to choose a name already 
in use); QED indicates the type of interaction. 

If there are new parameters related to coupling expressions, the user can edit a third file 
VariableName.dat. For example, to include the couplings between Z' and down quarks, 
the expression of the coupling appearing in the Feynman rule can be written as 



9z'dd = 



s w c w V 2 3 J 



I±^ + C 2 ^I^ (4.1) 
2 2 3c w 2 ' V ' 



where sw and cw are sin9w and cos 9\y, respectively. C\ and C<i are one for a SM like Z' . 
This can, of course, also be written as 

9z'dd ~ G i 2 + c 2 2 (4-2) 

where no link with the parameters in the SM is explicit. The declaration of these parameters 
appears as a list in VariableName.dat 

CI #first variable name 
C2 #other variable name 

The second step is to run the ConversionScript .pi. The script uses the above three 
files as input and creates the parameter card, param_card.dat and the file couplings. f, 
as well as other files needed by MadEvent. The third and final step is to include in these 
files the numerical values related to the model. The user model parameter card is in LHA 
format and is similar as the SM one, but also includes the information about the new 
particles and coupling parameters. The masses and widths of the new particles, as well as 
the new coupling parameters, need to be set by the user to their correct numerical values. 
In the couplings. f file, expressions for the coupling strengths have to be provided. The 
Perl script already takes care of the formats for the couplings in such a way that they are 
in HELAS compliant. For our example, the expression in Eq. ( [4.1[ ) can be implemented as 

GZPD(l)=dcmplx(Cl*( -ez*(-Half + sin2w/Three) ) , Zero) 
GZPD (2) =dcmplx (C2* ( -ey/Three) , Zero) 

The possibility of testing the couplings values and/or other parameters like the masses 
and the widths is also available. The program testprog can be compiled and run and it 
prints the values for all masses and couplings that will be used by MadGraph. A second 
program couplingsvalues writes out the names and the corresponding values of cou- 
plings in a format that can be read by external tools like the width and decay calculator 
BRIDGE |§. 
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5. Matching of jet production by parton showers and matrix elements 



For many processes at the LHC, the major backgrounds include multijet production, either 
pure QCD production or jet production in association with weak vector bosons or top 
quarks. Also when generating signal processes, such as Higgs boson production, it is often 
important to understand the jet activity, i.e. the probability for extra jet production, to 
be able to discriminate from the backgrounds. 

Parton showering, e.g. by Pythia, is well known to give a good description of parton 
emissions when the emitted partons are close in phase space. To describe hard and well sep- 
arated partons, matrix element generators such as MadGraph/MadEvent should be used. 
However, even in that case, it is necessary to perform parton showering and hadronization 
in order to get a realistic description of the event, in particular when detector simulation 
is needed. In order to combine these two descriptions, it is essential to use some kind of 
matching between them to avoid double counting of emissions in overlapping phase space 
regions. 

One such matching scheme is the CKKW algorithm [13, 14 1. This scheme uses a k± 
measure to separate emissions into two phase space regions, a low-A;_|_ region described by 
the parton showers and a high-&;_|_ region described by the matrix elements. In order to get 
a smooth transition between the two regions, the matrix element emissions are treated as 
similarly as possible to parton shower emissions. This is done by clustering the event using 
the k± jet clustering algorithm 47, 48] to find its corresponding "parton shower history", 
i.e., the sequence of parton emissions which would have been necessary for the parton 
shower to generate the event. The event is then reweighted by using the values of k± at 
the clustering nodes as scales in the running a s couplings (as is done in the parton shower), 
and a Sudakov suppression factor is applied to get the probability for this particular event 
to be generated without further emissions above the k± cutoff scale. After that the event 
is showered, allowing only shower emissions below the cutoff scale. 



An alternative way to accomplish an equivalent result has been proposed by Mangano 4£j. 
In his approach, the multiparton events are reweighted by a s factors in the same way as in 
the CKKW prescription, but no Sudakov reweighting is performed. Instead, the event is 
showered, and then discarded if the showering generates emissions harder than the phase 
space cutoff. 

In the current MadGraph/MadEvent release, both of these matching schemes are im- 
plemented. For the MLM matching option, the showering and vetoing of events with 
too large emissions are implemented in the Pythia interface included in the Pythia-PGS 
package (see Section |6.3|) , while the CKKW matching, including both a s and Sudakov 
reweighting, is so far implemented only at the matrix element creation level. 

In the original version of the MLM matching procedure, which is implemented in Alp- 
gen J7j , the phase space separation is defined in terms of cone jets. The MadGraph/MadEvent 
implementation allows the user to choose whether to use cone jets or for clustered jets, en- 
abling a more immediate comparison to the original CKKW formulation as implemented 
in Sherpa ||. 

A comparison between event generators implementing matching of matrix elements 
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and parton showers is underway fl9|| , including MadGraph/MadEvent, Alpgen, Sherpa, 
HELAC [5C] and Ariadne [51]. V^+jets production at Tevatron and the LHC is there 
used as a case study. 

It should also be noted that a similar matching method for matching of MadEvent 



matrix elements to Pythia showers has been implemented independently by Mrenna [15 
and has been successfully used to describe Tevatron VF+jets backgrounds. 

6. Tools 

The parton-level events generated with MadEvent are stored in the so-called "Les Houches 



accord event file format" [52]. Within this format enough information for each event is 
available so that some of the data analysis usually performed during the event generation 
phase, such as plotting, estimating PDF's errors or scale variations can be deferred to a 
later stage. To this aim, simple routines have been developed to perform some tasks "off- 
line", i.e., directly on the event files produced. The main reason is to improve versatility 
and save time. Generating events is a CPU expensive activity, which, in some cases, can 
take many hours. Therefore, it is not desirable to have to rerun codes only for making 
new plots or switching from one scale choice to another. Another important advantage in 
working directly with the events is that the tools developed are "independent" of how the 
events were generated and can be used with any event set in the Les Houches format. In 
this respect, the applications presented below can be used with events produced by any 
matrix element generator. The expert user is invited to develop his/her own tools and 
make them available to the physics community. 

6.1 Analysis platforms: MadAnalysis and ExRoot Analysis 

Two platforms for performing studies on the event samples at all stages of the simulations 
(parton, hadron and detector level) are available. Here we describe them briefly. 

MadAnalysis is a fortran-based analysis tool which allows you to select events (by 
setting further acceptance or selection cuts) and produce plots for events read in LHEF 
(parton level) and LHC Olympics (detector level) formats. Plots are simple ASCII files, 
output either in a minimal format (just title and comments and xy list) which can be 
read by GnuPlot |53| with a simple script (provided), or in TopDrawer [|54| format. A perl 
script to overlay plots of the corresponding quantities for two event sets (or the same set 
at different stages of the simulation chain) is also provided. 

Particles, objects and reconstructed objects (like jets) in the final state can be organized 
in classes, on which cuts and plots are based. Particles or objects in each class are ordered 
with respect to a one-particle quantity variable (the default is the transverse momentum). 
For instance, it is natural (at the parton level) to define any light quark and gluon as a jet. 
The user selects the cuts and the details of the plots by editing a ma_card.dat. Typical 
quantities that can be plotted are one-particle quantities, such as transverse momentum 
(p T ) and the rapidity (y) or pseudo-rapidity (77), or multi-particle quantities, such as two- 
body or three-body invariant masses (m(i, j),m(i, j, k)), or "distances" (such as R(i,j) or 
kT(i,j))) for each pair of final-state particles. A rather exhaustive library of kinematic 
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functions is available to the user, as well as the possibility of defining new functions of the 
momenta of any number of particles in a very user-friendly way. 

The ExRootAnalysis package allows to store and analyze events in a ROOT tree 
format [55]. Normally, ROOT allows to store data in several different formats. However, 



there are some crucial differences, which makes the ROOT tree format more attractive 
than others: 

• information is stored in arrays of objects, which enables efficient storage (compression 
using the ROOT gzip algorithm) and retrieval; 

• the same C++ classes are used for creating the ROOT tree and for analysing the 
stored data. 

ROOT tree objects are created from particles generated by MadGraph and Pythia or 
objects produced by PGS (in most of the cases physics objects like jets, electrons, etc.), in 
order to perform analysis in a ROOT environment. 

The ExRootAnalysis package can be subdivided into several subsystems: basic frame- 
work of few classes providing event loop, event selection and basic operations with a ROOT 
tree file; modules selecting events and objects to be analysed at per event and per object 
level; and modules analyzing selected events. 

For example, a selector module can select and group partons generated by MadGraph 
into several classes (such as leptons, jets, top quarks, etc.) according to their status and 
particle identification number. Any number of classes can be specified in the configuration 
file. After the partons have been classified, an analysis module can be used to produce 
series of standard plots for each class of partons. 

Documentation on the content of the ROOT tree is available on the web j56]| . 

6.2 Decaying unstable particles in the SM and beyond 

MadEvent includes a tool, Decay, that performs the decays of unstable particles in the 
Standard Model directly on the parton-level events (weighted or unweighted). At present, 
a total of 68 decay modes are included for r, W, Z, t and h decays. The advantage of using 
Decay is obvious: from the point of view of MG/ME the generation of events is faster for a 
simple final state. When a detailed knowledge of the spin-correlations is not needed, Decay, 
which is only keeping the diagonal terms in the spin-correlation matrix in the helicity basis, 
is an accurate and very efficient tool to get a multi-particle final state. A generalization and 



considerable improvement of Decay has been recently released by Meade and Reece [46], 
under the name of BRIDGE. BRIDGE can take an arbitrary model as input and compute 
the widths (and hence branching ratios) for all the kinematically allowed two-body and 
three-body tree level decays. It can also decay generated events, choosing decay modes 
randomly according to the branching ratios. 

6.3 Hadronization and detector simulation using Pythia and PGS 

One objective in the latest development of MadGraph/MadEvent has been to facilitate 
hadronization and detector simulation. We have therefore made a package available which 



-17- 



includes the hadronization package Pythia 6.4 [11] and the fast detector simulation package 
PGS 4 1 57], as well as the parton density package LHAPDF [j5lf| , a stripped version of 
CERNLIB, and the utility STDHEP [59] which is used for communication between Pythia 
and PGS. The package also includes main programs for running Pythia and PGS, and an 
interface for reading MadEvent files into Pythia. This package is available for download, 
or can be used directly in the on-line event generation to simulate fully hadronized and 
detector reconstructed events. 

The MadEvent-Pythia interface reads the Les Houches Event file output of MadEvent, 
and communicates model parameters such as particle masses to Pythia. The behaviour 
of Pythia is determined by the input file pythia_card.dat. In the case where it finds 
particles in the events which are considered massless in MadEvent but massive in Pythia, 
such as electrons or muons, the particles are given their Pythia mass, and momenta are 
redistributed to account for the changes. The interface also allows for matching of jet 
production by Pythia parton showers with multiparton samples generated by MadEvent, 
see Section |5[ 

The outputs of the Pythia main program are a binary file with the full hadronic event 
information in the STDHEP output format and a text file in the Les Houches Event file 
format with information on resonances, jets, leptons, and missing energy. Jet clustering 
for the latter file is performed either using a /cr-clustering routine or a cone jet routine. 

The PGS main program reads the STDHEP file from the Pythia generation, and 
performs detector simulation using the detector setup given in the input file pgs_card.dat. 
The output is a text file in the PGS4 LHC Olympics format with information on the trigger 
status of the event, jets, leptons and missing energy. 

All the output files can be read by the ExRootAnalysis package to get ROOT files 
and plots for the different stages. They can also be read by the Mad Analysis package for 
simple analysis purposes (see Section |6.1| ). The events are consistently numbered in the 
output files throughout the chain. 



7. Applications to hadron collider physics 

In this section we present several examples of the types of studies that can be performed 
with the MadGraph/MadEvent package. The aim is to show how easily the various new 
features can be used to perform signal and background analyses with both theoretical and 
experimental aims. Even though the discussion is kept concise for space reasons, we stress 
that all the results obtained in the following could stand in a dedicated publication and 
that some of them are original and presented here for the first time. 

7.1 Higgs search in pp — > h — ► W + W~: signal and backgrounds 

For a Higgs boson of moderate mass, 140 GeV < nth < 170 GeV, the golden discovery chan- 
nel is via gluon-gluon production and successive decay into a pair of leptonically decaying 
W bosons, Fig. ^. 
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Figure 2: The diagram correspond 
ing to the process pp — > h — 



This process has a large cross section and also 
special kinematics characteristics, notably the angu- 
lar distance between the leptons coming from the Ws 
decays, which can be exploited to improve the sig- 
nal over background ratio. Reducible backgrounds are 
dominated by processes which involve single top and 
tt production, which typically yield a larger number 
of jets in the final state and can be therefore sup- 
pressed by imposing a jet veto. Many studies have 
been performed on this channel and its backgrounds 
both at the experimental and the Monte Carlo level, 
sometimes with different generators, selection cuts and 

approximations 6Cj. 

In the following we discuss how such a study can be done using MG/ME. Our aim is to 
show how easy it is to generate all the necessary event samples in one go, from parton-level 
to detector objects, while keeping all the relevant non-trivial features of the signal (such 
as spin correlations) and the background (resonant and non-resonant contributions). In 
particular, the example makes use of the following features: 

• Inclusion of more than one process in the MadGraph generation (see Fig. ||). The 
different processes included in this way are then produced in proportion to their cross 
section, so the event file produced automatically gets the correct mix of unweighted 
events. 

• The possibility to generate the backgrounds corresponding to a given final signature 
(in this case e fjFbb), by including all the classes of diagrams (including resonant 
and non-resonant diagrams) consistently at once. The information on the relative 
contributions from the various resonances is available on event-by-event basis. 

• The interfaces to Pythia for decay and hadronization, PGS for detector simulation 
and MadAnalysis for plotting, analyzing and comparing observables at the various 
stages of the simulation. 



We consider that both VF's decay leptonically, and for simplicity we restrict to the 
different flavor case, e pF. As our purpose is illustrative we use (a simplified version 
of) the selection cuts used in Ref. fl6l| . we choose the angular distance between the two 
leptons in the transverse plane, A(j>, as the discriminating variable. As it is well-known, 
this variable encodes the fact that in the signal the two leptons tend to be produced close 
in phase space, due to the constraints coming from angular momentum conservation and 
the purely left-handed couplings of the W to the leptons. 



The selection cuts used in the A</> distributions at the parton and PGS level shown in 
Fig. [| are: 

• leptons: 25 GeV < p £ T < 50 GeV, \rf\ < 2.4, AR(£,£') > 0.4, and m(£,£') < 80 GeV; 
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# Begin PROCESS # This is TAG. Do not modify this line 



pp>h> e- mu+ ve~ vm 01 # Signal Process 

QCD=0 # max qcd order 

QED=4 # max qed order 

HIG=1 # max hig coupling order 

end_coup # the coupling list is over 

pp > e- mu+ ve~ vm/h (32 # Irriducible background pp > W+W-> e- mu+ ve~ vm 

QCD=0 # max qcd order 

QED=4 # max qed order 

end_coup # the coupling list is over 

pp > bb~e-mu+ve~vm/h @3 # Reducible background non-, single-, double-resonant 

QCD=2 # max qcd order 

QED=4 # max qed order 

end_coup # the coupling list is over 

done # the process list is over 

# End PROCESS # This is TAG. Do not modify this line 

# Begin MODEL # This is TAG. Do not modify this line 
heft 

# End MODEL # This is TAG. Do not modify this line 

# Begin MULTIPARTICLES # This is TAG. Do not modify this line 
p uu~dd~ss~cc~g 

# End MULTIPARTICLES # This is TAG. Do not modify this line 



Figure 3: The proc_card.dat used for signal — ► h — ► W + W~ — ► e~ fj, + 242/^ and the irreducible 
and reducible backgrounds. By a careful choice of the couplings and by vetoing the Higgs as an 
intermediate particle, the backgrounds are correctly included. The third process include non- 
resonant, single-top resonant and tt double-resonant production. 



• missing E T : E^' 1SS > 30 GeV; 

• jet veto: no jets with jP T > 30 GeV and \if\ < 3. 

We note that there is not much difference in the Acfi distributions in going from the 
parton-level events (where the Higgs has no transverse momentum and the two W's might 
have a very small one) to the fully simulated case. However, by plotting the pt of the 
leptons we checked that the initial state radiation has a non-negligible impact on the 
distributions, moving the average pt of the leptons considerably. 

Another very interesting feature to study is how the relative contributions coming 
from WbWb final states with top-quark non-resonant, single-resonant and double-resonant 
diagrams change as the jet veto is applied. This is a very good examples of a MC generation 



that has to be done carefully [62]. In fact single-resonant contributions are completely 
negligible and WbWb production is completely dominated by tt in generic areas of the 
phase space. On the other hand, the request of a jet veto dramatically enhances the non- 
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Figure 4: A<j> distributions of the two leptons in W + W events at the parton (left) and hadron 
(right) level. 

resonant and single-resonant contributions which account for more than 50% of the final 
WbWb background. 

7.2 Higgs CP properties: pp — ► hjj in the HEFT 

In this example we show how the Higgs effective theory implementation can be used to 
study the CP properties of a Higgs boson by looking at the angular distributions of the 
jets in the process pp — ► hjj |6^, |64}| . Our study is done at the parton level, but could be 
promoted to the hadron level without effort. 

We assume that the Higgs boson couples mainly to heavy quarks. In this case the 
Higgs boson will then be mainly produced by gluon fusion through a top quark loop. As 



discussed in Sec. [3^, for Higgs boson masses smaller than two times the top quark mass, 
rrift < 2mt, we can send the mass of the top quark in the loop to infinity rri; — > oo, to a very 
good approximation. Effectively, this means that we contract the top quark loops in the 
Feynman diagrams and get effective ggh, gggh and ggggh vertices. In Fig. || three diagrams 
with these effective vertices contributing to hjj production are presented. These vertices 
are implemented into the HEFT model of MadGraph, both for scalar and pseudo-scalar 
Higgs bosons. 

To investigate the CP properties of the Higgs boson we look at the angle A.(f>jj p5] , 
i.e., the angle between the transverse momenta of the two jets. We ask for one very forward 
and one very backward jet by applying the following cuts on the jets: 

p T (j) > 20 GeV, ARjj > 0.4, - rj h \ > 4, Vjl ■ rj h < 0. (7.1) 

These cuts lead to a signal comparable with Higgs boson production through W boson 
fusion. In Fig. ^, Acftjj is plotted with these cuts for a Higgs boson mass of Mh = 120 GeV, 
for the pure scalar and pure pseudo-scalar cases. 

7.3 Spin of a new resonance from lepton angular distributions 

Here we present the angular distributions for different intermediate particles in the process 
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Figure 5: Selection of the diagrams contributing to hjj production. The dashed line correpsonds 
to the auxiliary tensor particle T. For a pseudo-scalar Higgs the diagram on the rhs does not 
contribute. 
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Figure 6: A<j>jj distribution for pp — ► hjj production for a scalar and a pseudo-scalar Higgs boson 
through gluon fusion. The plots are normalized. Cuts are as in (0), M h = 120 GeV. 

pp — > X — > where X is a s-channel resonance. We consider the three cases where 

X is a spin-0, spin-1 or a spin-2 particle. 

This study exploits the possibilities of the user model framework to introduce new 
particles and interactions and also makes use of the full simulation chain, from parton level 
events to detector reconstruction. 

To minimize the effects of a non-zero transverse momentum of the intermediate state 
X, we use the angle 9 introduced by Collins and Soper to study spin correlations |66|] . The 
angle 9 is defined as follows. Let pa and pb the momenta of the incoming hadrons in the 
rest frame of the muon pair. If the transverse momentum of the muon pair is non-zero, 
then pa and pb are not collinear. Let us define an axis in such a way that it bisects the 
angle between pa and pb- The angle 6 is defined to be the angle between this axis and the 
fi + momentum in the muon pair rest frame. 

For (leading-order) parton level results, where we do not include extra radiation, the 
transverse momentum of the muon pair is zero. Then the angle 6 is the same as the more 
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commonly used angle 9*, i.e., the angle between the /x + momentum in the muon pair rest 
frame and the beam direction in the lab frame. 

In Fig. |7] the cosine of this angle is plotted for the three different cases of the spin of 
the particle X. This figure includes only qq initial states. 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



cos(theta) cos(theta) cos(theta) 

Figure 7: The normalized cross section as a function of cos£? in qq — > X — > Left for a 

spin-0 particle, center for a spin-1 particle, right for a spin-2 particle. No cuts applied. 

The spin-2 particle can also be created by gluon fusion, which dramatically impacts 
the distribution. The angular distributions of the muons depend, in general, on the way 
the X resonance was produced. In Fig. || the muons angular distribution is plotted for 
the gluon initial state and for the sum of the gluon and quark initial states, for a spin-2 
resonance with a mass of 1 TeV. 
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Figure 8: The normalized cross section as a function of cos 6* in gg — > spin-2 — > (left) and 

pp — > spin-2 — * )i + fj," (right). The mass of spin-2 particle is 1 TeV. No cuts applied. 

Including initial state radiation, showering and hadronization, does not modify the 
lepton distributions significantly. If we use PGS to simulate detector response we can get 
relatively reliable results for the CMS experiment at the LHC. In Fig. |9|, reconstructed 
events generated with PGS are plotted. In the same figure the parton level results are 
plotted with an acceptance cut for the rapidity of the muons |r/„+ „-| < 2.4. Due to this 
cut the events with 9 ~ (i.e., both muons in the same forward direction) and 9 ~ ir (i.e., 
muons in opposite forward-backward directions) are not detected. We see that for a 1 TeV 
resonance the lepton distributions hardly change at all. 
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Figure 9: The normalized cross section as a function of cos 9 in pp — > X — > for the 

reconstructed events after simulation of extra radiation using Pythia and detector response using 
PGS (dotted, red) and parton level events with rapidity cut \r)a+ a- I < 2.4 (solid, blue). Left for a 
spin-0 particle, center for a spin-1 particle, right for a spin-2 particle. Mass of the X particle is 1 
TeV. 
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Figure 10: The diagrams contributing to single-top and Higgs associated production. The first 
diagram is negligible in the SM and in type I 2HDMs while the second one is negligible in the 
MSSM and in type II 2HDMs. 



7.4 Single-top associated Higgs production in a generic 2HDM 

In this example we make use of the 2HDM implementation in MadGraph/MadEvent and 
of the associated calculator. 

The it associated Higgs boson production is well known to have a non-negligible cross 
section at LHC (~ 1 pb) for moderate values of the Higgs mass (m^ ~ 100 GeV). This can 
be estimated starting from the large it pair production cross section and then accounting 
for the large top-Higgs Yukawa coupling and phase space suppression. Since the single 
top production, in particular in the i-channel, also has a sizeable rate at LHC (more or 
less a third of the it production) it could naively be expected that the single-top asso- 
ciated production of the Higgs boson would be promising. It has been shown, however, 
that the cross section for single top in association with the Higgs is sizeably smaller than 
expected |)7[ |68| . 

The single-top associated production of the SM Higgs boson appears to be of the order 
of 100 fb for m/j ~ 100 GeV instead of the naive estimation of ~ 300 fb obtained by scaling 
the single-top cross section by a(tih)/a(ti). The main reason for this is a particularly 
strong destructive interference between the two dominant amplitudes, associated with the 
diagrams shown in Fig. [H]. It can be shown that each diagram contains a term proportional 
to m t which violates unitarity at high energies and cancels in their sum. 
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Figure 12: Cross sections of single top associated production of different Higgs bosons. On the 
left, the SM Higgs boson production cross section is shown together with the MSSM pseudoscalar 
^4° production cross section as a function of their masses for two different tan/3 values. On the 
right, the 2HDM pseudoscalar A production cross section is shown as a function of m H ± for type 
I and type II Yukawa couplings. A minimal pt of 20 GeV and a maximal rapidity of 2.5 is assumed 
for the jets. The factorisation and renormalisation scales are both set equal to m^. The b quark 
mass involved in the Yukawa coupling as well as its pole mass (in order to deal consistently with 
unitarity cancellations) are equal to the running mass at m^. The PDF used is CTEQ6L1. 




So, even though the sign of this interference can be 
guessed from unitarity requirements, its absolute magnitude 
at moderate energies, compared to the masses involved (e.g., 
the top mass), is surprising. The integrated squared ampli- 
tude for each of these diagram seperately is indeed up to 
three times larger than the total cross section. The total SM 
cross section at LHC as a function of is shown in Fig. 
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Figure 11: Extra diagram 

contributing to single top as Such a strong cancellation also takes place in the MSSM 

sociated Higgs production in 
the MSSM and 2HDM 111 case more as P ec ts have to be taken into account. 

To illustrate this, let us consider the decoupling regime where 
the mixing between h° and H° is small. For the light Higgs 
boson, the situation is similar to the SM. The WWh° coupling is close to the SM value 
and in the limit of large tan [3 where h° couples mainly to b, the bbhP coupling is close to 
the tth of the SM. The amplitude of the diagram where the Higgs boson is emitted from 
the initial state b quark appears to be slightly suppressed compare to the one associated 
to the diagram where it comes from the final state t quark (up to a factor 2 for identical 
Yukawa couplings), but the overall impact on the physical cross section stays small, or 
is even positive, since the negative interference previously described in the SM is also de- 
creased by the same factor. For the heavier scalar H° the coupling WWH° almost vanishes 
and one could then expect an enhancement of the total cross section. However there is an 
additional diagram involving a charged Higgs boson (see Fig. 11) which has to be taken 
into account due to the large W^H^H coupling. This diagram leads to a amplitude of 
roughly the same order and the same sign as the one involving only the W boson and 
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thus no particular enhancement is observed. The situation is similar for the pseudoscalar 
A (see Fig. 12) even though in this case the WW A coupling is strictly zero due to CP 



invar lance. 

Another possibility to get an increase of the cross section for the latter process would 
be to consider a 2HDM where the pseudoscalar A is relatively light and where the charged 
Higgs pair is much heavier (e.g., see |36|]) so that the negative interference cannot occur 



(i.e., the amplitudes associated with the SM like diagrams in Fig. 1C are dominant). This, 
of course, cannot occur in the MSSM where the masses of ± and A are linked at tree level 
through the mass relation m? H ± = m 2 A0 + m 2 v± . The resulting cross section for m^o = 40 



GeV is plotted in Fig. |12j both in case of type I and of type II 2HDM as a function of the 
charged Higgs mass. An enhancement at high m H ± is observed in the 2HDM type I case 
but the overall cross section stays much smaller than the SM one due to the reduced top 
quark Yukawa coupling. In type II models, the effect of varying mjj± is quite small. At 
low m H ± , the diagram involving the charged Higgs boson should contribute but its squared 
amplitude is more or less of the same order of magnitude as the negative interference it 
creates, so that its total contribution is negligible. Like in the MSSM, for all values of 
mjj±, the overall cross section appears to be slightly suppressed, more or less by a factor 
2, compared to the SM when the pseudoscalar is emitted from the initial state b quark. 

To conclude, the cancellations in single top associated Higgs production expected from 
the unitarity of the models considered together with the suppressions due to different 
emission configurations are very effective in reducing the cross section to a value close the 
SM one, which is probably too small to be successfuly measured at LHC. 

7.5 Comparison of strong SUSY pair production at the SPS points 

As a simple example showing some of the power of the SUSY implementation in Mad- 
Graph/MadEvent 4, we have chosen to compare results from inclusive strong SUSY pair 
production for the ten Snowmass (or SPS) benchmark parameter points [j69f |. In particular, 
the example makes use of the following new features: 



Implementation of the SUSY particle and interaction content [41] and the possibility 



to read SLHA files with mass, coupling and mixing information. 

• Inclusion of more than one process in the MadGraph generation. The different pro- 
cesses included in this way are then produced in proportion to their cross section, so 
the event file produced automatically gets the correct mix of unweighted events. 

• The possibility to define multi-particle labels also for scalar particles. 

• The interfaces to Pythia for decay and hadronization, PGS for detector simulation 
and ROOT for event analysis. 

7.5.1 Setup and generation 

To easily be able to see the relative importance of different groups of subprocesses, we 



generate the process using the proc_card . dat shown in Fig. 13. In this card we differentiate 
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between gluino pair production, squark of the first and second families pair production, 
third familiy squark pair production, associate production of gluinos and squarks, and 
associate prodution of third family squarks with first and second family squarks. In the 
generation all these subprocesses (in total 497 different subprocesses, since we make a 
distinction among the different squark flavors) are automatically generated in the correct 
proportions according to their relative cross sections. 

We next generate events at the LHC for each of the ten SPS points, and run the 
events through Pythia 6.409 [11], which decays the SUSY particles and performs parton 



showering and hadronization, and finally through PGS 4 [57]], a general-purpose detector 
simulation. We used the default LHC setup for PGS 4. The parameter cards for the points 



were generated using SoftSusy 2.0.5 |70|] and Sdecay 1.1a [Mf| and are available from the 



MG/ME 4 web sites. The results of the generation are: 

1. The total cross section for strong SUSY pair production. 

2. Cross sections for all the individual subprocesses, and summed cross sections for the 
different groups of subprocesses. 

3. Files with unweighted events at parton level, decayed and hadronized level and de- 
tector reconstruction level. 

4. ROOT files for event analysis at all three levels. 



7.5.2 Comparison of the results for the SPS points 

The cross sections and relative importance of the contributing subprocesses for the dif- 
ferent SPS points are presented in Table 4. The relative cross sections for the different 
subprocesses are calculated by MadEvent at leading order, but we should be able to take 
the results as a good indication of which processes are most important to take into account 
at each point. 

A striking feature of Table | is that between 68% (SPSla) and 92% (SPS2) of the 
total cross section is due to only the production of valence quark partners, gluinos, and in 
some cases top partners, constituting ~ 25 processes out of almost 500. This is of course 
due to the dominance of the valence u and d quarks over other quark species in the parton 
distributions of the proton, which is larger the heavier the squarks are. 

We also look at kinematical event distributions for the different points, after decay, 
hadronization and detector simulation. One of the most suggestive distributions is the 
total transverse event energy, Ht, defined as Ht = SPjJ* 8 + YJp^ ptons -|- since it gives 
an indication of the mass scale of the produced particles. In Fig. [14| the Ht distributions 
are shown for a reduced set of SPS points. We can clearly see how the position of the peak 
(or peaks) is correlated with the masses of the particles produced, see Table [|. For the 
ten SPS points, the peak positions are at about 65-85% of the sum of the masses of the 
produced particles, with a slightly increasing ratio as the masses increase. The exception 
is the very low mass stop pair production peak at SPS5, where the peak position is below 
40% of the sum of the produced stop masses. This is due to the large proportion of the 
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# Begin PROCESS # This is TAG. Do not modify this line 



pp > gogo @1 # Process 

QCD=2 # max qcd order 

QED=0 # max qed order 

end_coup # the coupling list is over 

pp > S1S1 @2 # Process 

QCD=2 # max qcd order 

QED=0 # max qed order 

end_coup # the coupling list is over 

pp > S2S2 @3 # Process 

QCD=2 # max qcd order 

QED=0 # max qed order 

end_coup # the coupling list is over 

pp > Slgo @4 # Process 

QCD=2 # max qcd order 

QED=0 # max qed order 

end_coup # the coupling list is over 

pp > S2go @5 # Process 

QCD=2 # max qcd order 

QED=0 # max qed order 

end_coup # the coupling list is over 

pp > S1S2 @6 # Process 

QCD=2 # max qcd order 

QED=0 # max qed order 

end_coup # the coupling list is over 

done # the process list is over 



# End PROCESS # This is TAG. Do not modify this line 

# Begin MODEL # This is TAG. Do not modify this line 
mssm 

# End MODEL # This is TAG. Do not modify this line 

# Begin MULTIPARTICLES # This is TAG. Do not modify this line 
p uu~dd~ss~cc~bb~g 

51 ulul~urur~dldl~drdr~slsl~srsr~clcl~crcr~ 

52 blbl~b2b2~tltl~t2t2~ 

# End MULTIPARTICLES # This is TAG. Do not modify this line 

Figure 13: The proc_card.dat used for the generation of all strong 2^2 SUSY processes. The 
processes are grouped such that gluino pair production is in the first group, first and second gen- 
eration squark pair production is in the second group, third generation in the third and associated 
production in the last three groups. Please note that further subdivisions have been made in 
Table § 

event energy carried away by the two LSPs, which largely balance in the detector so that 
the total missing energy detected is relatively small. If we add 2m x o to the peak positions, 
they reach between 85% and 100% of the mass of the produced particles, again with a 
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Table 4: Cross sections for the ten SPS points and the contributions, in %, from different sub- 
processes for each point. In the process definitions, q = ul.r, cIl,r- For gg and tt* production, the 
dominant channel is gluon fusion, with qq annihilation contributing to less than 10% of the numbers. 
The exceptions are point 2 and 5, respectively, where we have separated the gg and qq channels 
by the "+" sign. The last two columns show the contributions from associated gluino-(anti)squark 
production where the processes in column 4 are excluded, and (anti)squark-(anti)squark production 
where the processes in column 5, 7 and 8 are excluded. The NLO cross sections are calculated using 
Prospino2 0, §,§. 



H,. in strong SUSY pair production at the LHC 




Hp GeV/c 



Figure 14: Ht distribution for SPS points la, 2, 3, 5 and 9. The peak or peaks of the distribution 
is correlated with the masses of the produced SUSY particles. 



slightly rising trend going to higher masses. 

In the present study, we did not make use of the most important feature of Mad- 
Graph/MadEvent: its efficiency in calculating multiparticle final states. In particular, we 
used Pythia to perform the decay of the SUSY particles produced. For the scalar quarks 
this should be a reasonable approximation, since no spin correlations are expected (al- 
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SPS point 


la 


2 


3 


5 


9 


Peak position(s)(GeV) 


800 


1200 / 1900 


1300 


200 / 1100 


2200 


Particle masses (GeV) 


g: 600 
q: 550 
X?: 100 


g: 780 
q: 1550 
Xl : 120 


g: 930 
q: 830 
Xl : 160 


ii: 260 
3: 720 
g: 650 
X?: 120 


g: 1290 
g: 1260 
X?: 200 



Table 5: The positions of the peaks in the Ht distributions and the masses of the particles mainly 
produced, as well as the LSP, for the SPS points la, 2, 3, 5 and 9. 



though there might still be effects from interference or Breit-Wigner curve distortions in 
some parameter regions). However, if we want to study angular distributions where spin 
correlations might be important, the gluino decay should be done in MadEvent. For this 
kind of refinements, efficiency will start to become an important factor. It is therefore 
better to use different treatments for the different parameter points, in order to make an 
optimized choice for which processes to include at each point. So this should be seen rather 
as a preparatory study for such a more elaborate analysis. 



7.6 Inclusive W+jets matched samples: comparison with the Tevatron data 

At the Tevatron, as well as the LHC, many interesting signals include an isolated lepton, 
missing energy and hard jets. This means that weak boson production plus QCD jets 
is an important background. Until recently, this background was estimated using parton 
shower Monte Carlo's such as Pythia, where the largest jet multiplicity is + 1 jet, and 
additional jets are produced by parton showering. As discussed in Section || this gives a 
well-founded description of jets with small transverse momentum, and jets close in phase 
space, but for well separated hard jets it misses important non-logarithmic effects such as 
interference between diagrams. For that case, matrix element calculations are necessary. 
In order to simultaneously describe well-separated jets of different multiplicities, some kind 
of matching is needed between these descriptions. 

and Z boson production with multiple jets is for several reasons an excellent 
testing ground for matching procedures; there is a simple "central" 2 — ► 1 process, there 
is a hard scale set by the mass of the vector boson, and there is data from the Tevatron 
to compare to. A comparison between several implementations of matrix element-parton 
shower matching of W + jets production at hadron colliders is in progress fl9[| . 

In order to assess the robustness and flexibility of the method, there are several pos- 
sibilities for variations. Keeping to the principal matching method used in MadEvent and 
Pythia, a modified MLM method using the /c^-clustered jet definition, the most prominent 
variable is the so-called matching scale, used as cutoff between the matrix element and 
parton shower descriptions for jet production. Other natural variables are the scales for 
the running of a s , which can be varied in the parton shower and/or in the matrix element 
generation, and the scale used in the parton densities, which might also be varied indepen- 



dently. In Fig. 15, we have chosen a fixed matching scale of 10 GeV and varied the scale 
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Figure 15: The p± of the W boson and the parton-level differential 1 — > 2 jet rate for VK+jets 
production at the Tevatron, using matching of matrix elements and parton showers with MadEvent 
and Pythia. Please refer to the text for further discussion of the plots. 



for a s in the matrix element generation and the partons showers by a factor 2 up and down 
from the default scale, set by the k± in the clustering corresponding to that power of a s . 

Fig. 15 a shows the pj_ of the W boson in W production at the Tevatron, including 
hadronization of quarks and gluons but without underlying event simulation. The full black 
line shows the matched sum of the different jet multiplicity contributions, which are shown 
in different colors and line styles. For very low p±w the pure parton shower dominates, 
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but from the matching scale and up, the different jet multiplicities become more and more 
important, until around 150 GeV, where all the four included multiplicities are of similar 
order of magnitude. The band in the plot corresponds to adding a factor | — 2 in the 
arguments of a s in the matrix element generation only. It can be seen that the simulation 
with default parameters reproduces the shape of the Tevatron data j73|] (the normalization 
has been fitted by 20%). 

Fig. 15b shows the differential jet rate going from 1 — > 2 jets, after parton showers 
but before hadronization. It is included as an example, showing very clearly the transition 
between parton showers (in this case from the 0-jet and 1-jet samples), below the matching 
scale \0gQ2 = log A;™* = 1, and matrix elements (from the 2-jet sample and upward) above 
the matching scale, as well as the relative smoothness of the transition. The jet rate is 
the variable used to perform the matching, although here only particles with |?7| < 2.5 
are included in the jet definitions, explaining the low-end tails of the higher multiplicity 
samples. 



8. Conclusions and Outook 

In this paper we have presented the new release of the MadGraph/MadEvent package which 
aims at providing a tool to perform the simulation of both signal (Standard Model and 
beyond) and backgrounds within one framework. It is worth summarizing the key points 
of our approach. 

First, MG/ME is a user-driven Monte Carlo. Our package does not contain a library of 
processes but creates the code specific to the user's requests. It contains an easy-to- extend 
library of physics models which the user can choose from, including the Standard Model, 
the Higgs effective-theory, the MSSM, the most general 2HDM. 

Second, event generation with MadGraph is available from the web, in addition to the 
downloadable source code. All functionalities of our package can be accessed directly via 
a web interface without the need of any installation/compilation. User's code creation, 
event generation and detector simulation can be controlled by filling (or uploading) simple 
input cards and the results in the form of event files and plots can be downloaded from the 
user's personal database. MG/ME event generation is based on an algorithm which parcels 
the overall work in many small independent jobs. At present, our resources include three 
medium-sized PC clusters that are available and open to the public for code creation and 
the generation of (limited-size) event samples. The web server software is portable and can 
be easily installed on any cluster running PBS or Condor. Extension to work with the Grid 
software are under development. The source code for this, and support for experimental 
as well as theoretical groups is available upon request to the authors. 

MG/ME is designed to be a tool to ease the communication between theorists and ex- 
perimentalists. Any theorist, model builder or phenomenologist, can implement and care- 
fully test his/her own field theoretic model in MG/ME. The corresponding phenomenology 
at colliders, such as the LHC or an ILC, can be easily studied up to the simulated detector 
level by themselves and/or with the help of experimental colleagues. Both signal and back- 
ground can be simulated within the same framework. Eventually, the experimental groups 
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can pick it up and promote it to a "exp-grade" analysis by using their detector specific 
tools. 

The advantages of having "one framework for all" is manifest for both theorists (model 
builders) and experimentalists. For the former, it is important that the tool is very flexible, 
the learning curve is mild and no specific programming skills are needed. In particular, 
all the tedious and error-prone (and sometimes even difficult) tasks like matrix element 
creation, cross section integration and event generation, are automatically taken care of. 
On the experimental side, the insertion of a code in the simulation chain is normally a time- 
expensive process which involves testing and validation. For MG/ME this can be done, 
once for all, through a small set of simple Standard Model processes where comparison 
with available MC can be quickly done. 

The further developments that are planned for the package are also towards building 
a framework suitable for theorists/experimentalists interactions. 

In the current version, we have added a semi-automatic framework to implement new 
physics models. Even though this framework is quite general and very simple to use, it 
has two main drawbacks. First, the Feynman rules of a given model have to be available. 
This is trivial for very simple models but it can be become rather tedious and cumbersome 
for richer or more complete models (such as SUSY). Second, the range of models that can 
be easily implemented is still limited by the kinematic form of the interactions in HELAS, 
which has to be similar to those present in the Standard Model. We are currently working 
on a package that starting from a generic Lagrangian generates the Feynman diagrams in 
a MadGraph compatible format and on a more flexible version of HELAS. 

One feature that is under testing and will be soon included is the possibility of selecting 
the so-called decay chains in the generation of very rich multi-particle final states arising 
from the multi-staged decays of heavy particles in arbitrary models. By using a matrix 
element-based approach to decays, all information about possible correlations, such as those 
coming from the spin of the intermediate resonances, is available and can be exploited in 
the analysis. 

Another line of development is towards building general analysis tools that can max- 
imally exploit the information encoded in the matrix element to improve the accuracy or 
the sensitivity of measurements of key parameters (such as masses or couplings). These 
techniques have been successfully employed in several LEP and the Tevatron analyses, but 
always on a process specific basis. A general and automatic approach would be certainly 
welcome. 
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